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ABSTRACT 

The polarization signatures of the blazar emissions are known to be highly 
variable. In addition to small fluctuations of the polarization angle around a 
mean value, sometimes large (> 180°) polarization angle swings are observed. 
We suggest that such p henomena can be interpreted as arising from light-travel- 
time effects within an underlying axisymmetric emission region. We present 
the first simultaneous fitting of the multi-wavelength spectrum, variability and 
time-dependent polarization features of a correlated optical and gamma-ray flar¬ 
ing event of the prominent blazar 3C279, which was accompanied by a drastic 
change of its polarization signatures. This unprecedented combination of spec¬ 
tral, variability, and polarization information in a coherent physical model allows 
us to place stringent constraints on the particle acceleration and magnetic-field 
topology in the relativistic jet of a blazar, strongly favoring a scenario in which 
magnetic energy dissipation is the primary driver of the flare event. 
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1. Introduction 


Blazars are the most violent active galactic nuclei. They are known to emit nonthermal- 
dominated radiation from radio fr eque ncies to 7 -rays, with strong variability across the entire 
spectrum (e.g., Aharonian et ah 12007 ). Such emission likely originates from a relativistic jet 
that is directed close to our line of sight. The blazar spectral-energy-distribution (SED) has 
two components. It is generally acknowledged that the nonthermal radio through optical- 
UV emission is synchrotron radiation from ultrarelativistic electrons. In leptonic models, the 
high-energy component, from X-rays to 7 -rays, is usually interpreted as resulting from inverse 
Compton scattering by the same nonthermal electron population of soft seed photons. Seed 
photons can either come from the synchrotron component itself (synchrotron-self Compton 


= SSC; e.g., iMaraschi et al.l 11992c iMarscher fe Gear! 119851) . or from external photon fields 


Dernier et al. 
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). However, a hadronic ori 
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gin of the high-energy com- 

ponent cannot be ruled out (e.g 

, Mannheim & Biermann 

1992 


Miicke & Protheroe 

2001 ; 

Bottcher et al 

2013 

). Recently, 2 

jhang & Bottcher 

(2013f) have demonstrated that the high- 


energy polarization signatures can serve as a powerful diagnostic between the leptonic and 
hadronic processes. 

Despite intensive observational and theoretical efforts, the particle acceleration mecha¬ 
nism, the energy source for flaring activity, and the inner-jet physical conditions, such as the 
magnetic field topology, are not well understood. On the observational side, the synchrotron 
component is known to be polarized, with polarization degrees ranging from a few to tens of 
percent, in agree ment with a synchrotron origin in a partiall y ordered magnetic field. Sev¬ 
eral authors (e.g.. iLvutikov et al.ll2005t iPushkarev et al.ll2005l) have shown that the observed 
polarization degree and angle may reveal a helical magnetic field structure. Both the polar¬ 
ization degree and angle are known to be highly variable. Large (> 18 0°) polar i zatio n ang le 


swings h ave been freque ntly observed, and some recent observations (IMarscher et al. 2008 


2010 


Abdo et al . 2010 aj) have shown that they are sometimes accompanied by simultane¬ 


ous optical and 7 -ray flaring activities. When the polarization angle is not rotating, small 
fluctuations of the polarization signatures around a re latively stable mean value are often ob¬ 
served (e.g. lAleksic et al.ll20l4 iMorozova et al.ll2014h . Time-symmetric polarization profiles 
are observed in both cases, in particular a complete ~ 180° polarization angle swing is often 
accompanied by a drop of the polarization degree to nearly zero followed by a recovery back 
to its value before the PA rotation (see the above references). On the t heory side , the general 
formalism for synchrotron polarization is well understood (e.g.. iWestfoldl 1195911 . Neverthe¬ 
less, the spectrum, variability and polarization signatures have never been combined into one 
coherent physical model to understand the inner-jet physics. Models for synchrotron polar¬ 
ization necessarily take into account the magnetic field topology and polarization-dependent 
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synchrotron emissivity, but often apply simple, time-independent, power-l aw electro n popu - 
lations, and usually ignore the generation of high energy emission (e.g.. lLyutikov et all 20051). 
Thus they cannot produce full broadband spectral-energy-distributions (SED) and variabil¬ 
ity. On the other hand, models presenting detailed calculations of broadband emission and 
variability typically assume a chaotic magnetic held, and ignore any polarization-dependent 
emissivity and polarization information. Additionally, previous emission models hardly pro¬ 
vide any constrai nts o n the parti cle acceleration mechanisms and the energy source. In a 
recent paper, IZhang et ah (120141 ) (hereafter Paper I) have presented a general parameter 
study as a first step to combine the full SED, variability and polarization signatures in a sin¬ 
gle physical model. This paper has demonstrated that polarization-angle swings, correlated 
with high-energy haring activity, only require a dominant helical magnetic held structure in 
a straight jet, but no non-axisymmetric jet features, such as bends or helical how streamlines. 


In this paper, we will investigate additional geometric effects, namely the shape of the 
emission region and the magnetic held topology, on the synchrotron polarization signatures 
and multiwavelength variability. We find that the apparently time-symmetric polarization 
prohles can be obtained by an axisymmetric emission region along with full consideration of 
the light-travel-time effects. We will then apply the results obtained here and in Paper I to 
present the hrst simultaneous fitting of snapshot broadband SEDs, light curves in various 
frequency bands, and time-dependent polarization signatures, from a correlated optical and 
7 -ray haring event of the prominent blazar 3C279, which was accompanied by a drastic 
change of its polarization features. This unprecedented htting combination in one coherent 
model allows us to place stringent constraints on the particle acceleration and its energy 
source, as well as the magnetic held structure and its variation, in the inner jet of the blazar. 
Our results strongly favor a scenario in which magnetic energy dissipation is the primary 
driver of the hare event. We will describe our model setup in Section 2, present a general 
study of the geometric effects in Section 3, ht the 3C279 hare event in Section 4, and discuss 
the results in Section 5. 


2. Model Setup 

Our setup is designed to highlight the combined geometric and light-travel-time effects 
with the simplest physical assumptions. The underlying model considered here assumes a 
cylindrical emission region which travels on a straight trajectory defined by the jet boundary 
and encounters a hat stationary disturbance. The entire region is pervaded by a predomi¬ 
nantly helical magnetic held, with possibly an additional turbulent component. The origin 
of the disturbance could be a shock, which will convert bulk kinetic energy into random rno- 
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tion of relativistic particles to generate a flare; but it may also may also be associated with 
other mechanisms suc h as m agnetic reconnection, which will dissipate magnetic energy to 
accelerate particles (IGuo et al.l ( 20141 ). see Section 4). In the comoving frame of the emission 
region, the disturbance will travel through the emission region, and temporarily change the 
plasma conditions at its location,primarily the magnetic field structure and the nonthermal 
particle distribution through a modification of the acceleration time scale, and hence gen¬ 
erate a flare. As the disturbance leaves a given point in the emission region, the plasma 
conditions will recover back to the initial state. Most importantly, due to the light-travel- 
time effects (LTTEs), the observer will receive signatures of flaring activities induced by the 
disturbance, from points at different longitudinal distances along the jet at the same time. 
Owing to the axisymmetry of the emission region, this leads to time-symmetric polarization 
profiles (see Sections 3 and 4). 

The above model is simulated by the combination of the two- dimensional Monte-Carlo/Fokker- 
Planck radiation transfer (MCFP) code developed by Chen et ah ( 2012 . 20141 and the 3D 
multi-zone synchrotron polarization ray-tracing code developed by Zhang et ah] (12014 1. Fig. 

[Tj shows the model and code setup. The MCFP code considers a cylindrical emission region, 
assumes axisymmetry for all physical properties, and divides the region evenly into multiple 
zones in the radial and longitudinal directions. The electron distribution in each zone evolves 
in time according to a locally isotropic Fokker-Planck equation. The disturbance is also a 
cylindrical region, which will instantaneously modify the physical conditions of the zones at 
its current location, hence modifying the Fokker-Planck evolution. When it moves out of a 
zone, the physical conditions will instantaneously be restored to the initial state. Here, the 
physical meaning of an “instantaneous” change is that its time scale is less than one time 
step of the simulation. The code will then calculate the time-dependent emissivities based 
on the particle distribution at each time step, and apply a Monte-Carlo method to trace the 
photons and handle the Compton scattering. Therefore, all LTTEs are naturally included. 

All calculations are performed in the comoving frame of the emission region, and the resulting 
fluxes are properly transformed to the observer’s frame at the end of the simulation. 


The 3DPol code is focused on the synchrotron polarization signatures. It uses the 
same geometry and physical conditions, and the particle distributions as calculated by the 
MCFP code. However, since the polarization requires a 3D description, the emission region 
is further divided evenly into multiple zones in the </> direction, but the axisymmetry of all 
parameters is still kept. Since the environment of the emission region is generally considered 
as optically thin, and our treatment of the synchrotron polarization focuses on the high- 
frequency radio through optical/UV range, synchrotron-self absorption and Faraday rotation 
effects are neglected in this study. As in the MCFP code, all calculations are performed in 
the comoving frame. In each zone, we project the magnetic field onto the plane of sky in the 
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comoving frame, and use the time-dependent electron population generated by the MCFP 
code to evaluate the polarization degree ( PD ) at various frequencies. Since the net electric 
vector is perpendicular to the projected magnetic held direction, the electric vector position 
angle, also known as the polarization angle {PA), is easily obtained. In this way, we obtain 
the Stokes parameters (without normalization) at various frequencies at every time step in 
each zone via 

(/, Q, U) v = L„* (1, LQ cos 2 0 E , LQ sin 2 0 E ) (1) 

where L u and 1I„ are the spectral luminosity and the PD at frequency v, respectively, 
and 9e is the electric vector position angle for that zone. The code then employs ray¬ 
tracing to calculate the relative time delay to the observer for each zone, in order to take 
full account of the LTTEs. Since the emissions from different zones are incoherent, the 
total Stokes parameters can be obtained by linear addition of the Stokes parameters from 
each zone for emissions reaching the observer at the same time. Finally, as in the MCFP 
code, all emissions are properly Lorentz transformed (boosted) to obtain the luminosity and 
polarization information in the observer’s frame. 

Due to the relativistic aberration, even though we are observing blazars nearly along 
the jet in the observer’s frame (typically, 9* hs ~ 1/T, where T is the bulk Lorentz factor 
of the outflow along the jet), the angle 9 0 b s between LOS and the jet axis in the comoving 
frame it is likely much larger. Specifically, if 9* hs = 1/T, then 9 ohs = 90°. Paper I has shown 
that the polarization signatures are only weakly dependent on the exact value of 9 0 b s , thus 
in all of the following simulations, we choose 60x3 = 90° in the comoving frame, and hence, 
the Doppler factor <5 = (T [1 — /3 r cos6 l * bs ])~ 1 = T. 


3. Geometric Effects 


In this section, we use the quasar PKS 1510-089 as an example to investigate a num¬ 
ber of generic flaring scenarios as case studies. PKS 1510-089 is an Flat-Spectrum Radio 
Quasar (FSRQ). Its flaring activity in March 2 009 was well covered by multi-band monitor 


ing in the infrared, o pt ical, X-rays , and y -rays flAbdo et al.ll2010bl : iD’Ammando et al .112011 


Marscher et ah 2010). Chen et al.l (20121) presented model fits to snap-shot SEDs and light 


curves, in which they found that the most favorable scenario is an external Compton (EC) 
model with the external radiation held originating from a dusty torus. 


The purpose of this section is to study geometric effects on the polarization signatures, 
in particular the PA variability. The geometries are chosen in such a way that the first 
two cases correspond to 180° PA swings, while the other two cases are aimed to study to 
study small-scale PA fluctuations. In order to facilitate a direct comparison with paper I, we 
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Fig. 1.— Left: Sketch of the geometry used in the MCFP code, adapted from Chen et ah 


( 20121 ). Right: Sketch of the geometry (in the co-moving frame) of 3DPol. The model uses 
cylindrical coordinates, (r, </>, z), with n r , nn z being the number of zones in the respective 
directions. We define a corresponding Cartesian coordinate system (x, y , z) where z is along 
the axis of the cylinder and the x-axis is along the projection of the LOS onto the plane 
perpendicular to z. The Cartesian coordinates (x 0 ,y 0 ,z 0 ) are defined so that x 0 is along 
the LOS and zq is the projection of the cylindrical axis onto the plane of the sky. Both 
Cartesian coordinate systems are in the co-moving frame of the emission region. 0 O b s is the 
observing angle between x 0 and z. Consequently, if 9 ohs = 90°, ( x,y,z ) = (x 0 ,y 0 ,z 0 ). The 
cyan, dark-green and maroon regions represent far-side (left), middle and near-side (right) 
zones, respectively. 


choose the same quiescent state parameters, similar to those used in C hen et ah (120121 ). In 
addition, we keep the assumption that the disturbance is due to a shock. However, the shock 
parameters and the geometry may vary. Also, for the purpose of these generic case studies, 
we do not introduce any turbulent magnetic field contribution. Since both the MCFP and 
3DPol codes are time-dependent, we allow for an initial period for the electrons and the 
photons to reach equilibrium, before we introduce the parameter disturbance produced by 
the shock. In all our results the plots illustrate the results after this equilibrium has been 
reached. As the flaring activity exhibits different characteristics in duration and in strength 
for different cases, we define similar phases in the flare development for the purpose of a 
direct comparison. These phases correspond approximately to the pre-flare, early flare, flare 
peak, late flare, and post-flare states. 

We define the PA in our simulations as follows. PA = 0 corresponds to the electric 
field vector being parallel to the projection of the cylindrical axis on the plane of sky. An 
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increasing PA corresponds to counter-clockwise rotation with respect to the LOS, to 180° 
when it is anti-parallel to the projected cylindrical axis (which is equivalent to 0 due to 
the 180° ambiguity). Based on the above definition, the Stokes parameters normalized by 
luminosity for one zone with its projected magnetic field directed in the range of 0 - 45° 
is in the form of (1, — |g|, |w|), in the range 45° - 90° it is (1, |g|, |w|), for —45° - 0 it is 
(1, — \q\, — |w|), and for —90° - 45° it is (1, |g|, — |u|). This convention will be frequently used 
in the following text. All results are shown in the observer’s frame. 

Table 1 lists some key parameters. The emission region is a cylinder of length Z and 
radius R , traveling at a bulk Lorentz factor of T in the observer’s frame, when it encounters a 
flat, uniform, stationary disturbance of length L and radius A. The Fokker-Planck equation 
takes into account the acceleration time-scale f acc and escape time-scale t esc , both in units 
of Z/c, however, the total number of electrons is conserved, hence the escaping electrons 
will be balanced by thermal electrons picked up from the backgroimd. In the quiescent 
state the emission region is pervaded by a helical magnetic field Bh oriented at a pitch 
angle 9 b, and possibly a turbulent contribution, where the total magnetic-field strength is 
B. However, in this geometric effect study, we do not introduce any turbulence, therefore 
Bh = B. Initially a power-law distribution of nonthermal electrons with a power-index p 
and minimum and maximum energies 7 min and 7 max , respectively, with density n e , will have 
evolved to an equilibrium according to the Fokker-Planck equation in the emission region 
before the interaction with the disturbance. The disturbance will change the layers in the 
emission region at its location into an active state. In cases 1, 2 and 4, the disturbance, 
which is in the form of a shock, will compress the local magnetic field, so that the magnetic 
field strength will be amplified by a factor of B s /B, and the pitch angle will change to 0 S B . 
In case 3, the shock will instead shorten the acceleration time scale by a factor of tl cc /t acc , 
so as to increase the acceleration efficiency. 


3.1. B-Field Compression, Z < 2 R 

In this scenario, the shock will instantaneously increase the toroidal magnetic field 
component at its location, probably by compression, so that it will increase the total magnetic 
field strength and change its orientation in the affected zones. The new magnetic field will 
be kept until the shock moves out of that zone, when it instantaneously reverts back to its 
original (quiescent) strength and orientation, which can be attributed to dissipation. Paper 
I has considered a very extreme case, where the shock increases the toroidal component by 
a factor of ten, so that the total magnetic field is increased by approximately a factor of 
seven, and 6b increases from initially 45° to approximately 84°. Such strong alteration in 
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Parameters 

General Properties 

Bulk Lorentz factor F 

15.0 

Total (helical) magnetic field B ( G ) 

0.2 

Initial electron density n e (10 2 cm -3 ) 

7.37 

Initial electron minimum energy 7 m i n 

50 

Initial electron maximum energy 7max 

20000 

Initial electron spectral index p 

3.2 

Electron acceleration time-scale t ac c {Z/c) 

0.09 

Electron escape time-scale t esc {Z/c) 

0.015 

Orientation of LOS 0 o b s (°) 

90 


Parameters 

Case 1 

Case 2 

Case 3 

Case 4 

Helical pitch angle 6b (°) 

45 

45 

65 

45 

Length of the emission region Z (10 1G cm) 

8.0 

18.0 

18.0 

8.0 

Radius of the emission region R (10 lb cm) 

6.0 

4.0 

4.0 

6.0 

Length of the disturbance L (10 16 cm) 

0.8 

1.8 

1.8 

0.8 

Radius of the disturbance A (10 lb cm) 

6.0 

4.0 

4.0 

1.778 

B s /B 

\/3 + l 

\/3 + 1 

— 

\/3 + 1 

n 

75 

75 

— 

75 

^acc Aacc 

— 

— 

1/2.5 

— 


Table 1: Summary of model parameters. Top: General properties valid for all cases. Notice 
that 7 m i n , 7 m ax and p describe the injection of nonthermal particles and the electron distribu¬ 
tion will reach an equilibrium before interacting with the disturbance. Bottom: Parameters 
for each case. The s-superscript indicates the parameters during the presence of the shock. 
The parameter B s /B is the magnetic-field amplification factor, tlcc/^acc is the acceleration 
time scale shortening factor, and d s B is the magnetic-field pitch angle in the shocked region. 


the magnetic held results in a drastic change in the polarization signatures; however it will 
require a huge amount of energy conversion within small time windows when the disturbance 
is turned on and off, which is unlikely to happen in practice. Here we investigate a case with 
a moderate change in the magnetic held, where 9 b of the total magnetic held increases from 
45° to 75°. To mimic a spherical volume for the emission region, we choose R : Z = 1 : |, 
so that Z < 2 R, and the disturbance will occupy the entire layer, i.e., A = R (hereafter 
Geometry I). 

Since the total magnetic held is not increased much, we observe in Fig. [2] (upper left) that 
the synchrotron hux exhibits only a small hare. It is also obvious that the PD is photon- 
energy dependent (Fig. [2] lower left) since PKS 1510-089 has a visible external thermal 
component emerging at optical-UV frequencies, which diminishes the observed polarization 
percentage, especially in the UV. The dotted lines in Fig. [2] (lower left) show the intrinsic PD 
without the thermal contamination. The spectral curvature of the realistic electron spectrum 
from the detailed simulation based on the Fokker-Planck equation leads to an additional 
frequency dependence of the PD. The features above ~ 100 eV are due to the electron 
spectral cutoff. For the time dependence of the polarization signatures, the PD gradually 


























decreases to nearly zero at the flare peak, then recovers back to its initial state at the end 
of the flare. Moreover, the PA shows a continuous 180° swing, except for the radio/mm 
polarization, which is due to the extremely small flare amplitude at low frequencies. As is 
elaborated below and in Paper I, such polarization variability features can be explained as 
the combined effect of electron evolution and LTTEs. 

Initially, the directions of the projected magnetic field will cover the range of —45° to 45°. 
Due to the axisymmetry of the initial conditions, the initial \u\ component will be balanced 
out, hence the initial Stokes vector is (1,—|g|,0), leaving PA at 270° (equivalent to 90°). 
When the shock moves into the emission region, it will increase the toroidal component of the 
magnetic field at its location, resulting in stronger emission and enhanced synchrotron cooling 
(Fig. [3] left). Additionally, during the presence of the disturbance, the active region has 
6b = 75°. The combined effect is that the active region will have a larger contribution to the 
polarization, with the Stokes parameters generally in the form of (1, |g|, |w|) or (1, |g|, — |w|). 
Due to LTTEs, although the shock is flat in the co-moving frame, the active region at equal 
photon-arrival times will be distorted, as shown in Fig. b. In the early flare phase, 
only the region facing the observer (red in Fig. [4]) contributes to the flare, where the initial 
Stokes parameters in the form of (1, — \q\, |w|) will be replaced by the active state (1, \q\, |w|). 
Therefore the initial PA = 270° will gradually move to 225°, and the PD will decrease at 
the same time. When the flare gradually rises up to its peak, emission from the far side 
of the cylinder can be seen, which possess negative u component. This will diminish the 
positive u from the right side of the cylinder, meanwhile the initial negative q has mostly 
be canceled out. Consequently the PD will continue to drop, and the PA will rotate from 
225° to 180°. When the flare reaches its maximum (green in Fig. |4 Jd) , for a short period the 
shock-enhanced fluxes from both sides of the cylinder will be comparable. As a result, the u 
component will be canceled out, leaving a positive q in the active region and a negative q in 
the quiescent region. If the shock is adequately strong, the polarization contributed by the 
active region will be higher, leading to a positive q for this short period. Therefore, the PD 
drops to nearly zero for all bands, while the PA reaches 180° for infrared, optical and UV. 
However, the PA rotates back to the initial 90° for the radio, as its flare amplitude is not 
strong enough to dominate over the quiescent emission. After the peak, the flaring region 
moves to the far side, so that the polarization signatures gradually revert back to the initial 
state in time-symmetric patterns. 
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Fig. 2.— Case 1: A moderate change of magnetic field in Geometry I. Upper left: Snap-shot 
synchrotron SEDs, including the external photon field. SEDs are chosen at approximately 
the pre-flare (black solid), early flare (red dashed), peak (purple short dashed), late flare (blue 
dash-dotted) and back to equilibrium at the post-flare (orange dash-dot-dotted) states, with 
the dotted line for the external photon field contribution. All SEDs are chosen with the same 
time bin size. Lower left: PD vs. photon energy for the same time bins as the SEDs in the 
top panel, with the external photon contamination considered, where dotted lines represent 
the PD without the contamination. Upper right: PD vs. time with external contamination, 
at radio (black solid), infrared (red dashed), optical V band (purple short dashed), UV (blue 
dash-dotted). Lower right: PA vs. time for the energy bands as in the top panel. 

3.2. B-Field Compression, Z > 2 R 

In order to isolate the dependence of the polarization signatures on the geometry of 
the emission region, we studied a scenario in which we kept all the parameters exactly the 
same as in the previous case, except for the size of the emission region. We choose an 
emission region with the same volume, but now R : Z = | : 3 (hereafter Geometry II). The 
general electron spectrum is similar to the previous case, except minor differences due to the 
electron cooling rate, escape rate, etc. introduced by the different geometry. However, we 
notice that the flare amplitude is slightly increased, thus the radio polarization signatures 
now behave similarly to the other bands. The major difference in this case is that in this 
scenario, the PD goes through two minima near zero with an elevated plateau inbetween. 
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Fig. 3.— Electron spectra at different time steps for a moderate change in the magnetic field 
(left) and a shortening of the acceleration time scale (right). 0: Pre-flare equilibrium; 2: cen¬ 
tral shock position; 4: shock just leaves the emission region; 8(11): post-shock equilibrium. 
The regions between the magenta vertical lines represent the electron energies that corre¬ 
spond to the photon energies we choose in the polarization vs. time plots: Solid corresponds 
to radio, dashed to infrared, short-dashed to optical, and dashed-dotted to UV. 

The PA exhibits a similar step-like pattern with a plateau at 180°, although, in total, it 
still completes a full 180° swing. The reason is that unlike in the previous case, during the 
peak, the entire equal-photon-arrival-time ellipse is maintained within the emission region 
for a finite amount of time, during which the polarization signatures remain unchanged. 


3.3. Acceleration Efficiency, Z > 2 R 

We investigate one additional scenario where the shock leads to more efficient particle 
acceleration, by instantaneously shortening the acceleration time scale at its location. In 
this case, we choose the initial toroidal component slightly dominant, with the application of 
Geometry II (we have shown one case with Geometry I in Paper I). Due to the more efficient 
acceleration, the electrons are accelerated to higher energy, and the spectrum becomes harder 
(Fig. [3] right). Therefore, the synchrotron component exhibits a considerable flare at higher 
energies. However, at lower energies the electron spectrum stays almost unchanged. Since 
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the magnetic field remains unchanged, the radio band shows no time dependency of either 
flux or polarization signatures. 

At the beginning of the flare, the luminosity of the active region is drastically enhanced; 
as a result, the observed PD will mostly be attributed to the red region in Fig. [4] (right). 
In this small region, the magnetic field is well ordered, consequently we observe that PD 
can shoot up to over 60%; meanwhile PA quickly rotates to about 210°, which represent the 
6b in this region. Similar to the previous case, during the flare peak, the entire emission 
region becomes axisymmetric, therefore PA rotates back, and both PD and PA exhibit a 
step phase. After the peak, PD and PA revert back to the initial state in time-symmetric 
patterns. 


3.4. B-Field Compression, Localized shock 

We present one case for the combined geometric effect in Fig. [71 This time we apply a 
moderate change in the magnetic orientation and strength during the presence of the shock, 
with the same volume of the emission region as in Geometry I (R : Z — 1 : |), but now the 
shock will not occupy the whole layer of the emission region, but only the center of it, so 
that the region that will be affected by the shock, has A : Z — | : 3, (hereafter Geometry 
III). The current situation is equivalent to a Geometry II case but surrounded by a large 
quiescent region. Hence the positive q in the shock region will not be able to dominate over 
the negative q in the quiescent region. Consequently, we find a similar behavior as in Case 
3, where the PA rotates back to 90° and exhibits a step phase. The difference is that the 
PD slightly decreases during the flare peak, since the shock will give a boost to the toroidal 
contribution to the polarization, diminishing part of the initial poloidal contribution. 


4. Application to 3C279 


Based on the insights gathered in the above study and Paper I, we now present the first 
simultaneous fitting of snapshot SEDs, multifrequency light curves, and time-dependent 
polarization signatures of a blazar, using the example of a flare of 3C279. The FSRQ 3C279, 
located at a moderate redshift of z = 0.536, is one of the most well-observed members of the 
blazar class. It caught the attention of the high-energy astrophysics community due to its 
very bright gamma-ray flaring at the beginning of the Compton Gamma-Ray Observatory 
mission in the early 1990s, and h as since then been the target of many dedicated multi¬ 


wavelength observing campaigns flAbdo et al.l l2010al: lHartman et al.l Il996l; iMaraschi et al. 
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1994 Wehrle et al.lll998lh It is bright and variable on a large range of time-scales, across 
the electromagnetic spectrum, from radio through gamma-rays, and is one of only three 
FSRQs detected in very-high-e nergy gamma-rays (i.e., photon energies of E > lOOGeV) by 
the MAGIC telescope ( Albert et al. 2008 ). 


During an active (i.e., high-flux) phase from Nov. 2008 to Mar. 2009, 3C279 exhibited 
multi- wav elength flaring activity, with a period of optical polarization variations lasting ~ 20 
days f Abdo et al. 2010ah . The multi-wavelength light curve data are displayed in Fig. |8l 
They suggest that this flaring episode is actually composed of two sequential flares. The 
first flare dominates around MJD 54880, where a sudden doubling in gamma-ray flux is 
accompanied by a relatively small infrared-to-optical flare and rather erratic changes of the 
PD and PA. The second flare, which dominates the last ~ 15 days, constitutes a correlated 
flare of the infrared-to-optical and gamma-ray emissions, during which the PD drops to 
nearly zero and then recovers to ~ 20%, accompanied by a ~ 180° rotation of the PA. This 
was the first time that such a clear correlation between optical/gamma-ray flaring and a PA 
rotation was observed. Based on the apparently time-symmetric profile of the second flare, 
we suggest that it actually lasted ~ 20 days, but is overwhelmed by the decaying phase of the 
first flare for the first ~ 5 days. Here we present a consistent interpretation of the spectral 
energy distribution (SED), multi-wavelength light curves, and time-dependent polarization 
features, for the second flare. 

Table E] lists the most relevant parameters in our simulation that yields the best fit to 
the data of 3C279, and the results are shown, in comparison with the data, in Fig. El The 
low-frequency component of the SED is dominated by synchrotron emission from nonthermal 
electrons, along with a weak thermal component from the central accretion disk around the 
black hole powering the relativistic jet. We assume a leptonic origin for the high-energy 
component of the SED, which is composed of an SSC contribution, dominating from X-rays 
to soft gamma-rays, and an EC contribution which dominates the emission in the Fermi 
(~ GeV) range. We also assume that there exists a turbulent magnetic field component, 
which is initially uniform everywhere. We assume the disturbance to be associated with a 
magnetic energy dissipation process, e.g., magnetic reconnection, instead of a shock. It has 
been suggested that magnetic reconnection can dissipate a large fraction of the magnetic 


ener gy to produce a nonthermal power-law distribution of relativistic particles flGuo et ah 
2014I ). A possible underlying physical picture is that on the trajectory of the emission 
region, it encounters a flat stationary region, where the poloidal component of the helical 
magnetic field is in the opposite direction of that in the emission region. Therefore the 
poloidal component will be dissipated during the presence of the disturbance, due to magnetic 
reconnection, leading to particle acceleration (we use a simple injection to mimic this effect) 


and a locally stronger turbulent magnetic-field component (IDanghton et al.ll2011t iGuo et al. 
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2014). 


The time-dependence of the polarization signatures is similar to Case 1 of the previous 
section. Before interacting with the disturbance, the entire emission region contributes 
uniformly, and therefore, due to the axisymmetry of the underlying geometry and the mildly 
stronger poloidal component {9 b = 33°) in the quiescent state, the polarization is dominated 
by the poloidal contribution. When the interaction with the disturbance starts (red in Fig. 
ED, we assume that the poloidal component in this area decreases, but the toroidal component 
remains nearly unaffected, hence the quiescent-state polarization will be gradually canceled 
out by the flaring toroidal contribution. As a result, the PD drops and the PA starts 
to rotate towards the dominant toroidal direction. As the disturbance moves forward, the 
flaring region becomes larger, resulting in increasing fluxes in both the synchrotron and the 
Compton emissions. Towards the end of the interaction period, the flaring region (green 
in Fig. ED reaches its maximum size, leading to the observed flux maximum. Meanwhile, 
the magnetic field throughout the green flaring region in Fig. El has approximately equal 
contributions from the near and far sides of the cylinder, mimicking the initial axisymmetry, 
except for a mildly stronger toroidal magnetic-field contribution {9b = 62°). This flaring 
toroidal contribution is just strong enough to dominate over the poloidal contribution in 
the quiescent state. Consequently, the PD drops to nearly zero, while the PA reflects the 
dominant toroidal magnetic-field direction. After this flux peak, the flaring region gradually 
becomes smaller and moves to the upper-left (blue in Fig. ED, and the light curves and PD 
recover to their initial states with approximately time-symmetric profiles. However, since 
the toroidal component on the far side of the cylinder is opposite to that on the near side, 
the PA instead completes a ~ 180° swing to a direction that is indistinguishable from its 
initial position due to the 180° ambiguity of the PA. 


Discussions 


Larionov et al. 

2013; 

Morozova et al. 

2014) 


the observed > 180° PA swings are accompanied by one or several sequential apparently sym¬ 
metric PD patterns which drops from an initial value to zero then reverts back. In addition, 
both the PD and PA patterns appear to be smooth. Several mechanisms have been proposed 
to interpret the PA rotations, such as an emission region moving along a curved trajectory 
or a bending jet, or streamlines followin g helical magnetic-f i eld lines, or stochastic activation 


of in dividual zones in a turbulent jet flAbdo et al.l l2010at iMarscher et al.l |2010[ iMarscher 


2 014 ). While they may have their own virtue in understanding polarization variations, none 
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of these models has so far been able to explain spectral variability properties, symmetric 
light-curve profiles, and correlated symmetric polarization variability features in one coher¬ 
ent model. On the other hand, the LTTEs coupled with an axisymmetric emission region 
in our model naturally explain the apparently time-symmetric features of multi-wavelength 
light curves and polarization variations and their intrinsic correlations, which also appears 
to be the simplest model with the smallest number of fine-tunable parameters. 


Based on our model, the continuous 180° PA swings will pose some serious constraints 
on the physical background of the emission region, namely: the pitch angles of the quiescent 
and the active states shall be on each side of a critical point 9 C , a moderate strength of 
flare, a moderate change of the pitch angle of the helical magnetic held, and a geometry of 
the emission region where Z > R. If the pitch angles of the quiescent state and the active 


state are both within the range of 0 to 9 C or 6 C to 90°, where in our case 6 C ~ 55°, then 
the projected magnetic held will always provide the total Stokes parameters in the form 
of l,—\q\,u or l,\q\,u, respectively, resulting in the PA fluctuating around 270° or 180°, 
respectively. Additionally, if the hare amplitude or the pitch angle change is too weak, then 
throughout the haring activity the total polarization signatures will be dominated by the 
quiescent state, therefore PA will not complete a 180° swing. While if on the other hand they 
are too strong, then the active region will dominate the polarization signatures for a much 
longer time, giving rise to a step-like feature as shown in Paper I. Notice that, however, this 
feature is different from the “step-phase” feature introduced by a Z > 2 R geometry (Case 
2), because the step-like feature caused by strong variations can only reach 180° at the very 
peak of the hare, where the u component of the Stokes parameters can be fully canceled 
out. On the other hand, the “step -p hase” due to the geometry can be diagnosed in the 
observation as completely hat fllkeiiri et al.l 120111). 


In addition to the above constraints, a simultaneous fitting of all spectral, light curve, 
and polarization properties can exclude a wide range of possible scenarios. Therefore, our 
fitting results place unprecedented constraints, most of which do not depend on the details 
of the model. Since in a leptonic model interpretation, the gamma-rays from 3C279 are 
produced by the EC process, an excess of either external photon held or nonthermal electrons 
is necessary to produce a gamma-ray hare. However, the former fails to produce the data, 
since such excess is unlikely to generate a correlated, time-symmetric multi-wavelength hare, 
but instead may lead to an anti-correlated behavior in the synchrotron emission due to excess 
radiative cooling. Also, the infrared-to-optical hare amplitude is smaller than that of the 
gamma-ray hare, and the X-ray emission, which represents the low-energy end of the SSC 
emission, shows almost no variability. This implies that the total magnetic held strength has 
to decrease during the hare. We suggest that the dissipated magnetic energy can be converted 
into nonthermal particle energy through magnetic reconnection. Most importantly, as we 
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have discussed above, a smooth PA rotation and the nearly zero PD at the flare peak (green 
in Fig. [4| imply that the toroidal magnetic-field component should be mildly stronger than 
the poloidal one in the active state. However, since the excess of nonthermal electrons and 
the amount of the magnetic field alteration, which are linked by the reconnection process, 
have been well constrained, the poloidal and toroidal B-field components should remain 
comparable throughout the process. Moreover, the relative percentage of the turbulent 
magnetic-field contribution increases in the active state. Nevertheless, we notice that, just 
like most multiwavelength models of blazar emission, our model underestimates the radio 
flux, and also overestimates the PD there. This is because the radio emission is known to be 
produced on larger scales, on which the magnetic field appears more disordered. In addition, 
the environment in which our emission region is located is still opaque to radio emission. 


Models of relativistic sh ocks propagating through the jet have been widely used to ex¬ 
plain blazar flaring activ ities (IBottcher fe Dermerll201Q[ iJoshi fe Bottcherll2007l; iMarscher fe Gear 
19851 : ISpada . et al. 200 11). Such models naturally provide excess nonthermal electrons. How¬ 
ever, they are expected to compress the magnetic field, leading to a change in 9 b with an 
increase in the overall magnetic-field strength. Consequently, shocks are unlikely to fit our 
constraints. Alternatively, the fitting constraints strongly favor a magnetic energy dissipa¬ 
tion process during the flare. We find in our simulation that the dissipated magnetic energy 
during the disturbance is comparable to the necessary amount of particle energy increase re¬ 
quired to generates the flare. Simulations of magnetic energy dissipation have demonstrated 
that the energy stored in magnetic sh ear can be efficiently converted into a power-law distri¬ 
bution of relativistic particles flGuo et al.l 2014). This process will reduce the magnetic field 
component that is subject to dissipation and can therefore change the magnetic-field pitch 
angle. Moreover, the magnetic-field topology inside the dissipation zone is likel y to beco me 


turbulent, thus it will strengthen the turbulent magnetic-field contribution ( Daughton et al. 

201 lh . 


For the smaller fluctuations in the polarization signatures (generally PA varies by 
< 90°), although some a pparently symmetri c profiles can be no ticed, the general patterns 
appear very complicated (. Torstad et alJ 2013; Covino et ah 2015 ). In addition, those fluctu¬ 
ations often happen during the quiescent states, with lower PD. This impli es that so me in¬ 
homogeneity and more complex geometry, or some turbulence are required (IMarsc her 2014). 


We suggest that, if a step phase along with time-symmetric variability of the polarization 
features is observed, especially in the PA profile, as it will be less affected by a possible 
turbulent B-field component, this may serve as a powerful constraint on the size and the 
geometry of the emission region. Since the step phase will serve as a measure of the length 
of Z — 2 R, while the two symmetric PA (PD) alterations will imply the length of A. In 
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this way, the size and the geometry of the region that is affected by the disturbance can be 
constrained, while the size of the entire emission region can then be estimated by the flux 
and polarization percentage. 
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Fig. 4.— a) Sketch of the interaction between the emission region and the disturbance 
in the comoving frame of the emission region, at different epochs. The emission region is 
pervaded by a helical magnetic field and a turbulent component (only the helical component 
is sketched). The disturbance is stationary in the observer’s frame, but in the comoving frame 
of the emission region, the disturbance is then moving up with Lorentz factor T. The part 
of the disturbance that first encounters the emission region (shaded ellipses, the “front”) is 
the location of the injection of relativistic particles. The orange, red, green and blue regions 
refer to the locations of the disturbance before the flare (to), the rising phase (ii), peak (t 2 ) 
and declining phase (£ 3 ), respectively, b) The red, green and blue shapes indicate the shape 
and location of the flaring region, corresponding to the disturbance at t\ ~ t 3 , respectively, 
observed simultaneously, taking into account the LTTEs. c) Sketch of the projection of 
the helical magnetic field onto the plane of sky in the comoving frame. The upper panel 
illustrates the quiescent state, the lower panel the active state. The cyan, dark-green and 
maroon arrows represent the left side, center and the right side of the emission region shown 
in a), corresponding to the color regions in Fig. [T] respectively. || (_L) denotes components 
parallel (perpendicular) to the bulk motion direction. The two dashed lines indicate ±45°. 
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Fig. 5.— Case 2: A moderate change in magnetic field in Geometry II. Panels and line styles 
are as in Fig. [2] 



Fig. 6 .— Case 3: Shortened acceleration timescale in Geometry II. Panels and line styles 
are as in Fig. [2j 
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Fig. 7.— Case 4: Moderate change of the magnetic held with Geometry III. Panels and line 
styles are as in Fig. [2 
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Fig. 8.— Data and model fits to multi-wavelength SEDs, light cur ves and po la rizati on 
signatures durin g the s econd flare of 3C279 in 2009. Data are from lAbdo et al .1 fl2010al h 
Havashida et ahl (120121 ). Hollow (filled) d ata points refer t o the first (second) flare, a) 


SED, black squares are from Period E in [Havashida et al . (2012) (MJD 54897 - 54900), 
corresponding to the end of the second flare. The black curve is the model SED from the 
simulation at the same period; the red curve is the simulated SED at the peak of the flare. 
Hollow magenta circles are from Period D (MJD 54880 - 54885), corresponding to the end 
of the first flare. b,c) J and R Band flux: black squares are from the period of the second 
flare, hollow magenta circles are from the first flare; the black curves are the simulated light 
curves, d) Gamma-ray photon flux light curve: black squares and hollow magenta circles 
are 3-day averaged data from the second and the first flare, respectively, while green squares 
and hollow blue circles are the 7-day averaged data, respectively, e, f) PD and PA vs. time: 
black squares and hollow magenta circles are from the second and the first flare, respectively; 
curves are the simulated polarization signatures. 
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Parameters 

General Properties 

Bulk Lorentz factor F 

20 

Length of the emission region Z (10 17 cm) 

4.8 

Radius of the emission region R (10 17 cm) 

2.4 

Length of the disturbance L (10 17 cm) 

1.28 

Orientation of LOS 0 o bs (°) 

90 

Electron acceleration time-scale t acc (Z/c) 

0.36 

Electron escape time-scale t esc ( Z/c ) 

0.062 

Electron background temperature T (m e c 2 /ks) 

100 

Parameters 

Quiescent 

Active 

Total magnetic field strength B (10 —2 G) 

7.2 

5.4 

Helical magnetic strength, Bh (10 — ‘ 2 G) 

4.8 

2.6 

Helical pitch angle Oq (°) 

33 

62 

Initial electron density n e ( cm~ 

5.5 


Injected electron minimum gamma 7 m in,inj 


2000 

Injected electron maximum gamma 7 ma x,inj 


4000 

Injected electron power-law index pj n j 


4 

Injection rate Qi n j (10 43 erg * s —1 ) 


5 


Table 2: Summary of parameters. The parameters are defined in the same way as in Section 
3, except for a few differences. In this fitting, we assume the radius of the disturbance is 
equal to that of the emission region, so that A = R. As we mentioned in Section 3, the initial 
electron distribution will evolve to an equilibrium according to the Fokker-Planck equation, 
before interacting with the disturbance, therefore we only list the initial electron density n e , 
as the spectrum will be determined by the Fokker-Planck equation. We assume that the 
escaped electrons will be balanced by the electrons picked up from the thermal background, 
with a temperature T. The magnetic field has two components, a helical component Bh 
along with a turbulence; the total magnetic held strength will be B. The disturbance will 
change the layers in the emission region at its location into an active state. In such situation, 
the magnetic held energy will be dissipated, so that both the total strength and the helical 
strength will decrease and the pitch angle will alter. The dissipated energy will become the 
energy resource for particle acceleration, which we handle it by an injection at the front of 
the disturbance with an energy injection rate Qi n j into the emission region, with minimum 
and maximum Lorentz factors r y r nin,inj an d lmax,inji and power-law index p in j. Both the 
initial and the injected electrons will evolve according to the Fokker-Planck equation. 
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